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We suggest a method of generating distillable entanglement form mixed states unitarily, by uti- 
lizing the flexibility of dimension od occupied Hilbert space. We present a model of a thermal spin 
state entering a beam splitter generating entanglement. It is the truncation of the state that allows 
for entanglement generation. The output entanglement is investigated for different temperatures 
and it is found that more randomness - in the form of higher temperature - is better for this set up. 
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I. INTRODUCTION 

In quantum information we use quantum mechanics 
to perform information tasks in ways exceeding the ca- 
pabilities of classical systems 0. A natural and impor- 
tant question is then - what is it that gives us the power 
in quantum information and where does it come from? 
It has been shown by Knill and Laflamme that it is 
possible to make computations with exponential speed 
up over classical algorithms for certain tasks, using only 
one pure qubit and a resource of entirely mixed states. 
Bose et al. 0] have proposed a scheme where only one 
pure state is needed to generate entanglement between 
an atom in a pure state and a thermal field - no mat- 
ter what the temperature. Concerning mixed states and 
nonlocality, in Q Filip et al. used mixed states to violate 
Bell's inequalities. 

In optics, for example, we know that in the infinite 
dimensional case of the harmonic oscillator, if a state can 
be described as a statistical mixture of Glauber states, 
this state cannot be used to generate entanglement using 
a beam splitter |^. A true maximally mixed state 
cannot be used to generate entanglement via any unitary 
transformation whatsoever. This is obvious since if the 
input to an unitary is proportional to identity, so must 
its output be, in which case the state remains unchanged. 

Here we wish to address in part the question of how 
much mixedness we can deal with and still get entan- 
glement, and how we can do this. We use the simple 
idea of change in the dimension of occupied space to give 
an example of how we can get around even maximally 
mixed inputs to generate entanglement. By expanding 
the occupied Hilbert space, the overall entropy remains 
the same, and the state is no longer a maximally mixed 
state for that space, and thus may have some entangle- 
ment. Now, such an operation is not unitary, in the sense 
of the input Hilbert space. However, if the input state 
is viewed truncated state of a higher-dimensional 
space on which the unitary acts, and the unitary in ques- 
tion can expand the occupied Hilbert space, we can get 



unitary generation of entanglement. 

To illustrate this we use the model of finite trun- 
cations of a thermal state of optical modes entering a 
beam splitter. Such a state, truncated to dimension d, 
is exactly analagous to a thermal spin state, with spin 
S — {d — l)/2. We refer to these as thermal spin states. 
Indeed, in principle, we are not restricted to optics, and 
can imagine using real spin systems (please see Sec. IVTl . 
The beam splitter is the example we use of a unitary 
transformation which increases occupied Hilbert space. 
It is exactly this, along with truncation, that allows 
the generation of entanglement. In the case of infinite 
temperatures this gives unitary entanglement generation 
from a truncated maximally mixed state. 

We note that one must be careful with considering 
physical applications of this model. It is difficult to imag- 
ine a beam splitter type transformation acting on a finite- 
spin system, and even then, on the restricted space of 
the input, this is not a unitary transformation. In some 
cases infinite-dimensional systems can be modelled by 
considering finite spaces, but this is not what we are in- 
terested in here, since it is the truncation of the space 
that leads to the stark difference between the finite and 
infinite cases. We present this model simply as in illus- 
tration of the ideas using a real unitary that expands the 
occupied Hilbert space. The problems and possibilities 
of physical implementations of this model will be looked 
at in Sec. 

We begin in the first section by describing the beam 
splitter transformation and defining the measure for en- 
tanglement we use. Next we look at the entanglement 
generation for a truncated maximally mixed state. We 
then look at thermal states, again in the truncated sense. 
Finally we show a simple, if inefficient, protocol proving 
that the entanglement is distillable for all states consid- 
ered, then end with discussion and conclusions. 



II. BACKGROUND 
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In optics, the beam splitter is defined by its unitary 
action on a two mode Fock state, |m, n), returning the 
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output state 0, 



1+n) 



Ubs\m,n)= ^ f{m,n,M)\M,m + n- M), (1) 

M=0 



where 



f{m,n,M) = 



min(m,,M) 



E 



n\m\y/M\{m + n - M)\ 



p\{p - M + n)\{m - p)\{M - p)\ 



and where T and R are the complex transition and reflec- 
tion coefficients respectively, with normalisation |Tp + 
= 1. We use this as our definition of the beam split- 
ter on the general number state (i.e. not restricted to 
optics). 

We notice that for any input state whose density ma- 
trix is diagonal in the two mode Fock basis, the output 
entanglement is not effected by the phases of R and T, 
since it can be absorbed into a phase on the basis states. 
In all our states this is the case, so we do not need to con- 
sider the phase, and assume, without loss of generality 
that R and T are real. 

As a measure of entanglement we take the logarithmic 
negativity Ej^r, defined for a given state p as. 



E^{p) = log2||p^-||i 



(3) 



where pi are the eigenvalues of p"^-* , the partial transpose 
of p in subspacc A. For a state on a Hilbert space H = 
Ti-A ® 'Hb 7 the partial transpose in subspacc A is defined 
by its matrix elements 



(m, jslp^-^lfcylj^B) = {kA,jB\p\iAjB), 



(4) 



with d = 25* -I- 1 . This corresponds to a thermal state at 
infinite temperature. The output is then, 

Pout = UbsPinUl^ 

2S 2S (m+n) (jn+n) 

= d^EE E E f{^.n,M)f{m,n,M') 
x\M,m + n-M){M',m + n-M'\. (6) 



For any such input state of finite dimension we get an 
entangled output. Thus from "truncated" maximally 
mixed input states we have unitarily derived entangle- 
ment. The reason this is allowed is that the unitary op- 
eration does not act only on the space of the input; in this 
sense the beam splitter is not unitary on the truncated 
Hilbert space of the input. 

Figure n is a plot of entanglement against 5* for R = 
T = 1/ ^/(2) for truncated maximally mixed input states. 
We see an initial rise followed by a slow tail. We know 
that in the infinite S limit the entanglement must fall to 
zero, since this is a statistical mixture of Glauber states, 
which docs not entangle through a beam splitter (such 
states arc also called "classical" states) 
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where IiaJb) = \i)A ® \j)B & 'Ha ® Ti-B form an or- 
thonormal product basis. The logarithmic negativity is 
a widely used measure of entanglement and is an en- 
tanglement monotone ^3 ■ We note that a positive value 
of this measure does not necessarily imply the existence 
of useful, i.e. distillable entanglement (however, for our 
cases we prove distillability by independent means). 



III. MAXIMALLY MIXED (THERMAL) STATES 

We are now ready to consider unentangled, maximally 
mixed states, input to our beam splitter. 
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FIG. 1: Entanglement of two non entangled, maximally mixed 
states passing through beam splitter against S for R = T = 
I/V2. Here we see an initial peak followed by drop. In the 
large 5* limit this goes to zero - since this is nothing but a 
"classical" state of the harmonic oscillator. 

We may consider the initial rise in entanglement sur- 
prising since in some sense the higher the 5*, the closer 
we are to a "classical" state, and so we might think we 
should get lower entanglement. The rise can be explained 
however; as we increase the size of the spin of the input, 
the available space for entanglement increases - the max- 
imal entanglement follows ln{d), where d = dimension. 
Then, as more higher levels get occupied, they begin to 
destructively interfere and the entanglement begins to 
fall with S. Similar trends were observed in |0, [l3l| . 
where the entanglement was studied for spin coherent 
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states sent through a beam sphttcr and thermal states in 
the one-dimensional Heisenberg Model respectively. 

We can also consider optimality of entanglement gen- 
eration in terms of the beam splitter reflectivity R (since 
we are deling only with real values, T is set by R through 
|rp = 1 — Figure 121 shows the entanglement gener- 

ated by maximally mixed input states against S and 
For low S there are two maxima, at = 1/2 they are 
around 0.2. As S increases the peaks get closer and for S 
larger than around 3 they merge into one at i? = l/\/2. 
Note, however, that even if we chose R for each S to give 
the maximum entanglement, we will still see an initial 
peak in entanglement as we go from S" = 1/2 up. 
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FIG. 2: Entanglement of two maximally mixed states passing 
through two arms of a beam splitter against S and We 
see that the maximal depends on S and is at |_Rp = 1/2 for 
large S. The same trend is seen for all different temperatures 
of inputs - as long as their are two input beams, i.e. neither 
are vacuum (if this is the case, there is only one maxima at 
= 1/2 for all S). The lines are present simply to help 
show the texture of the plot. 



IV. THERMAL STATES OF DIFFERENT 
TEMPERATURES 



We now consider a general thermal state of our trun- 
cated system. We assume the Hamiltonian 



IS 



H = nfiw\n) (yij. 



(7) 



n=0 



giving the equalibrium thermal density matrix 



where Z = ^^^^ e 



ihoj/KT\ 



n){n\ 



(8) 



-nhuj/KT 



is the partition function. 



For two such thermal states incident on a beam splitter 
at different temperatures. 



FIG. 3: Entanglement of various thermal states against S of 
the input states for R = T = 1/^/2. The top four, unfilled, 
data sets correspond to when one port has a thermal state 
entering it and the other only the vacuum. The bottom four, 
filled, data sets are when both ports have thermal states en- 
tering them. 
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where Zi = J2 
beam splitter transformation is then, 



the output state after the 



Pout — UbsPinUl^ 

, 2S 2S 

Z1Z2 ^ ^ 

m— n— 

(m+n) (m+n) 

^ ^ f{m,n,M)f{m,n,M') 

M=0 M'=0 

x\M,m + n- AI) {Al\ m + n- M'\. (10) 

Figure |31 shows some plots for various different input 
temperatures, both for the case of a thermal state enter- 
ing one port and the other remaining empty (i.e. vacuum 
state, we can also think of this as zero temperature), and 
for thermal states entering both ports. We see that the 
general trends are the same for all temperatures, that 
of an initial peak followed by a slow decline. All plots 
are for reflectivity R = T = I/V2. This is chosen since, 
when plots similar to figure |21 were done for the different 
temperature inputs, the same trends were found, thus 
the overall trends are not affected. 

In general the case where the vacuum enters one port 
gives more entanglement - this can be explained as re- 
sulting from the fact the two beams do not interfere with 
one another, which destroys entanglement. 

Starting from the maximally mixed state, which rep- 
resents an infinite temperature thermal state, the height 
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of the peak is less for lower temperatures and the peak 
occurs earlier. We might expect the entanglement to be 
greater for lower temperatures, and that the maximally 
mixed state gives the lowest entanglement, since these 
are the most disordered and classical like states. The ob- 
served trend can be explained by noticing that for lower 
temperatures the higher dimensional states arc not as 
populated, restricting the possible entanglement. For the 
zero temperature we have the ground state which offers 
no entanglement. We also see a sharper fall after the 
peak for lower temperatures. 



V. DISCUSSION 

For all the output states mentioned here we can de- 
vise a protocol to distill a minimal amount of entan- 
glement. From equations 1)10(1 and © we can sec that 
if local projections arc made on both arms, onto the 
subspace spanned by the states |0)(0| and |4S')(4S'|, the 
remaining state is an entangled state of the form p = 
F|0, 0)(0, 0| + (1 - F)(|0, 45) + |45, 0))((0, 4S\ + (45, 0|), 
which has entanglement for any finite F. For finite S, F 
is also finite (hence we have entanglement). This can be 
considered as a two dimensional state and it can be easily 
shown that it has negative partial transpose, which, for a 
two dimensional, bipartite state implies distillablc entan- 
glement ^3 . As S goes to infinity, F goes to one and no 
entanglement can be found in this way as we expect. All 
measurements where the projection falls onto the remain- 
ing space are discarded. This scheme shows that for any 
finite S we indeed have distillable entanglement for all 
the thermal states, though it is not necessarily inefficient 
and it destroys at least some, of not most entanglement. 

Normally one might assume that from a maximally 
mixed state, l/(i(8)l/c?, no entanglement can be extracted 
via unitary operations. We see here that this is not the 
case and that when considering finite dimensions, there is 
an extra subtlety allowed by expanding of contracting the 
occupied Hilbert space. In some sense we can consider 
the operation as expanding the "mixcdness" of the in- 
puts into two, individually less mixed outputs with some 
entanglement. More precisely we have the input output 
entropy relation S (ptni) + S {pin2) < S{pouti) + S{pout2), 
where the inequality comes from the entropy that has 
gone to create entanglement (the overall entropy is con- 
served since the operation is unitary). From this per- 
spective it is not so surprising that more entanglement 
is generated for the more mixed (higher temperature) 
thermal input states. Conversely wc can think of the in- 
verse operation as contracting all the mixcdness into a 
separable mixed state of lower dimension. This must be 
true of many unitaries other than the beam splitter, in- 
deed we may consider the power of a unitary in terms 
of expanding or contracting the occupied Hilbert spaces. 
For example wc can imagine easily a unitary in a six di- 
mensional space acting on maximally mixed state of two 
qubits to create entanglement. The beam splitter simply 



represents a simple, well known and physical example of 
this typ e of unitary operation. The scheme by Browne et 
al. [13 , represents an application of a beam splitter used 
to extend the accessed Hilbert space to distill gaussian 
states. 

The morel we use here, as mentioned in the introduc- 
tion, is seemingly not readily very implementable. This 
model is intended as an illustration of the ideas using a 
real physical unitary that expands the occupied Hilbert 
space, rather than a proposal for an experiment. How- 
ever, it is hoped that that it could be useful as a basis 
for more sophisticated models. 

For this, one can think of the input being either a trun- 
cation of an infinite space, or a real thermal state of a 
finite system. With infinite dimension, for these ther- 
mal states to be physical, such a truncation should be 
physically imposed - this could happen for example, for 
a finite number of photons in a mode, or atoms in a 
Bose-Einstein condensate (BEC). A BEC gas with at- 
tractive interactions, for example, has a natural limit on 
the number of particles, imposed by the interactio ns |l6l . 
The beam splitter interaction also exists for BEC ■ A 
finite maximally mixed state for example, could be en- 
visaged as one of a pair of finite atom-number BECs in a 
maximally entangled state. In principle, then, this could 
offer an implementation of the scheme. 

One naive possibility may be the projection of a real 
optical thermal state onto a finite-dimensional subspace 
and then passing through a beam splitter, giving entan- 
glement. In such a case, the problem of entanglement 
generation becomes one of truncation or projection. This 
is related to the idea of measurement-induced entangle- 
ment (see e.g. ^^^3)- these schemes a measurement 
is made at the end of the protocol which gives guaran- 
teed entanglement conditioned on certain outcomes - it is 
the measurement that generates entanglement (whereas 
we would measure first and generate entanglement after). 
These schemes too implement projection, which does re- 
duce the size of the occupied space. The difference is 
that such schemes are concerned with guaranteed entan- 
glement generation, rather than exploiting the dimension 
of the occupied Hilbert space to deal with mixed state in- 
puts. 

Another possibility might be to map a mixed state 
of some finite dimension to a higher-dimensional space, 
and perform a unitary analogous to the beam splitter in 
expanding the occupied Hilbert space, creating entangle- 
ment. This could be, for example, a spin system in a 
thermal state mapped to the photon number states in 
a cavity, then passed through a beam splitter. Feasible 
mapping of spin systems onto light has been studied re- 
cently (e.g. [lalSfl)- In this sense we could begin with 
a genuine thermal spin state (even infinite temperature) 
and extract entanglement via unitary transformation. 
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VI. CONCLUSION 

111 this paper we have given an example of distiUable 
entanglement generation from a maximally mixed state 
via unitary operation. This result is explained by the fact 
that the state is in fact a truncated maximally entangled 
state, or equivalently, that the operation of the beam 
splitter is not unitary in the finite dimensional Hilbcrt 
space of the input. We have then looked at the entangle- 
ment generation for different thermal state inputs. Again 
against intuition, we found that the higher the input tem- 
perature, the more the resulting entanglement. This was 
put down to the fact that for lower temperature the popu- 
lation of the higher dimensional states is lower and hence 
the possibility for entanglement generation reduced. 

The source of the entanglement, or "quantumness" for 
this scenario can either be thought of as the truncation of 
the input state or the power of the unitary to expand the 



occupied Hilbert space. There are many natural ques- 
tions which arise, such as the minimal unitary way to 
extract entanglement from a maximally mixed state (as 
mentioned, the beam splitter is probably not the optimal 
unitary) . More comparisons of these finite systems to op- 
tics should lead to more understanding and intuition. As 
experiments get closer to high dimensional, finite sys- 
tems, (e.g. 12 if), these questions become more and more 
valid and useful. 



Acknowledgments 

This work was sponsored by Engineering and Phys- 
ical Sciences Research Council, the European Commu- 
nity, the Elsag spa, the Hewlett-Packard company and 
the Asahi Glass Foundation (Natural Science Research 
Assistance). 



[1] M. A. Nielsen and I. L. Chuang. Quantum Computation 
and Quantum Information. Cambridge University Press, 
2000. 

[2] E. Knill and R. Laflamme. Phys. Rev. Lett, 81(5672), 
1998. 

[3] S. Bose, I. Fuentes-Guridi, P. L. Knight, and V. Vedral. 

Phys. Rev. Lett, 87(050401), 2001. 
[4] R. Filip, M. Dusek, J. Fiurasek, and L. Mista. 

Phys. Rev. A, 65(043802), 2002. 
[5] M. S. Kim, W. Son, V. Busek, and P. L. Knight. 

Phys. Rev. A, 65(032323), 2002. 
[6] Wang Xiang-bin. Phys. Rev. A, 66(024303), 2002. 
[7] R. A. Campos, B. E. A. Saleh, and M. C. Teich. 

Phys. Rev. A, 40(1371), 1989. 
[8] G. Vidal and R. F. Werner. Phys. Rev. A, 65(032314), 

2002. 

[9] J. Lee, M. S. Kim, Y. J. Park, and S. Lee. J. Mod. Opt, 
47(2151), 2000. 
[10] J. Eisert. Ph.D Thesis, University of Potsdam, 2001. 
[11] J. Eisert, S. Scheel, and M. B. Plenio. Phys. Rev. Lett., 



89(137903), 2002. 
[12] D. Markham and V. Vedral. Phys. Rev. A, 67(042113), 
2003. 

[13] M.C. Arnesen, S. Bose, and V. Vedral. Phys. Rev. Lett., 

87(017901), 2001. 
[14] M. Horodecki, P. Horodecki, and R. Horodecki. 

Phys. Rev. Lett, 78(574), 1997. 
[15] D. E. Browne, J. Eisert, S. Scheel, and M. B. Plenio. 

Phys. Rev. A, 67(062320), 2003. 
[16] C. A. Sackett and R. G. Hulet. J. Opt B, 3(R1-R28), 

2001. 

[17] Y. Castin and J. Dalibard. Phys. Rev. A, 55(4330), 1997. 

[18] M. D. Lukin. Rev. Mod. Phys., 75(457), 2003. 

[19] A. S0rensen and K. M0lmer. quant-ph/0304008. 

[20] K. Hammerer, M. M. Wolf, E. S. Polzik, and J. L Cirac. 
quant-ph/0409109. 

[21] N. Gisin, L Marcikic, H. de Riedmatten, W. Tittel, and 
H. Zbinden. 2003. Proceedings of the 6th international 
conference on Quantum Communication, Measurement 
and Computing (QCMC 02), quant-ph/0301181. 



